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A theoretical method for treating collisions in the presence of multiple potentials is developed 
by employing the Schwinger variational principle. The current treatment agrees with the local 
(regularized) frame transformation theory and extends its capabilities. Specifically, the Schwinger 
variational approach gives results without the divergences that need to be regularized in other 
methods. Furthermore, it provides a framework to identify the origin of these singularities and 
possibly improve the local frame transformation. We have used the method to obtain the scattering 
parameters for different confining potentials symmetric in x,y. The method is also used to treat 
photodetachment processes in the presence of various confining potentials, thereby highlighting 
effects of the infinitely many closed channels. Two general features predicted are the vanishing of 
the total photoabsorption probability at every channel threshold and the occurrence of resonances 
below the channel thresholds for negative scattering lengths. In addition, the case of negative ion 
photodetachment in the presence of uniform magnetic fields is also considered where unique features 
emerge at large scattering lengths. 

PACS numbers: 32.80.Gc, 34.50.-s, 34.10.+X 


I. INTRODUCTION 

Frame transformation theory constitutes an important 
theoretical toolkit describing Hamiltonian systems with 
several potential terms which are dominant in different 
regions. A particular class of frame transformations, i.e. 
the local frame transformation (LFT), can provide piv¬ 
otal insights for quantum mechanical systems which are 
approximately separable, usually in different coordinates, 
at small or large distances but are not separable over all 
space. In their seminal papers, Fano[I| and Harmin[2-Q 
developed the principles of the LFT theory permitting 
the theoretical treatment of highly excited alkali atoms 
in the presence of external electric fields and interpreting 
on physical grounds the corresponding photoabsorption 
spectra. 

In addition, the unique properties and intuitive in¬ 
sights make the LFT theory the key framework in the 
investigations of atom-electron collisions under the ac¬ 
tion of uniform magnetic Q or electric fields The 

description of the corresponding photoabsorption spec¬ 
tra in terms of the LFT approach is accurate even though 
only a few channels were involved in the actual calcula¬ 
tions. 

However, in 1966 it was shown that the inclusion of all 
closed channels yields non-trivial effects when a charged 
particle is scattered by a zero-range potential in the 
presence of a magnetic field 0- In particular, it was 
shown that there always exists a bound state for neg¬ 
ative scattering length that is collectively supported by 
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all the closed channels (Landau levels) induced by the 
magnetic field. The properties of this bound state were 
subsequently investigated in the context of negative-ion 
photodetachment in weak magnetic fields m- However, 
similar effects also appear in ultracold two-body colli¬ 
sions in tight waveguide geometries yielding confinement- 
induced resonances of the Fano-Feshbach type0. Such 
scenarios have been studied in terms of LFT theory 
either for atom-atom scattering [l3l - [0 or dipole-dipole 
collisions [T^. LFT theory provided the most complete 
theoretical framework in order to describe the relevant 
collisional physics encompassing the concept of infinitely 
many closed channels. Interestingly, this permitted the 
theory to go beyond previous studies where the im¬ 
pact of the closed channels is included in an averag¬ 
ing sense [il,0. However, the integration of the con¬ 
cept of infinitely many closed channels in the LFT ap¬ 
proach yielded divergences, which required regulariza¬ 
tion tec hniq ues to be employed to remove the singular 
behavior [la Il4. [la [l^ . Evidently, an approach that in¬ 
corporates the concept of infinitely many closed channels 
while avoiding the infinities of the LFT treatment could 
be useful. Such a formalism could provide quantitative 
insights into the origins of the LFT singularities while 
preserving the physical intuition of LFT. Also, it could 
allow the systematic study of non-trivial effects arising 
from the full Hilbert space in collisional systems, e.g. 
the scattering of an electron in a Rydberg state from a 
neutral perturber[2^. 

In this _ paper, a Schwinger variational iF-matrix 

approach |2l|-l23| is developed for treating Hamiltonians 
with two potential terms which dominate at different 
length scales. We restrict our treatment in this paper 
to the situation of symmetric confining potentials in x,y 
and no force in z; the more general case of asymmet¬ 
ric potentials and/or forces in z will be treated else- 



2 


where. Note that the Schwinger variational principle 
has been successfully used for the quantitative study 
of Rydberg molecules (^. Rydberg atoms [1^, electron 
scattering [2^ [23 and electron-molecule collisions [28l - l^ . 
However, in most of these studies, the trial function is 
composed from a superposition of basis functions (often 
Gaussian orbitals), which requires numerical evaluation 
of complicated volume integrals. In our analysis, such 
complications are avoided due to the length scale sepa¬ 
ration between the two potential terms, permitting us to 
treat the scattering information analytically. 

This paper will focus on the case of a short-range spher¬ 
ically symmetric interaction potential with a confining 
potential that bounds all the degrees of freedom apart 
from one at large distances. Note that these types of 
Hamiltonians refer to the physical systems of atom-atom 
collisions in tight waveguides or negative ion photode¬ 
tachment processes in external fields. Under these con¬ 
siderations, the calculations within the Schwinger vari¬ 
ational framework are compared to the corresponding 
regularized LFT derivations for two different confining 
potentials (a harmonic oscillator and an infinite square 
well potential), yielding identical results. However, in the 
current Schwinger variational iF-matrix approach, the in¬ 
finitely many closed channels associated with the confin¬ 
ing potential produces finite results without regulariza¬ 
tion. 

Our analysis includes several applications of the cur¬ 
rent formalism to collisions in the presence of either har¬ 
monic or infinite square well confining potentials where 
the corresponding photoabsorption spectra for s-wave 
photoelectrons possesses unique features. Namely, the 
total photoabsorption probability vanishes at every chan¬ 
nel threshold regardless of the specific details of the con¬ 
fining geometries. This threshold phenomenon illustrates 
the impact of the closed channel physics on the colli- 
sional complex. In addition to the cases of harmonic 
and infinite square well confining potentials, the current 
development is applied also to infinite rectangular po¬ 
tentials or off-centered scattering in an infinite square 
well. These particular systems illustrate the physics of 
infinitely many closed channels in the absence of degen¬ 
eracies in the spectra of the confining Hamiltonians. As 
a final application, negative ion photodetachment in the 
presence of a uniform magnetic field is treated. Note 
this particular system has been thoroughly investigated 
both experimentally [3l|, [s^ and theoretically [1, [I^ [s^ . 
However, in most of the theoretical investigations only 
the open channel physics was considered and this permits 
us to explicitly show the impact of the closed channels. 
More specifically, our results are compared to the approx¬ 
imate LFT treatment in Ref. @ where only the open 
channels physics were considered. Qualitative changes 
occur only near thresholds for small scattering lengths. 
However, these features are enhanced the cases where the 
scattering length is comparable to the confinement length 
scale which implies enormous scattering lengths for lab¬ 
oratory strength magnetic fields or ultrastrong magnetic 


fields for the usual scattering lengths. 

Section H introduces the basic concepts which are used 
in the rest of this work. Section HI focuses on the deriva¬ 
tion of the Schwinger variational iF-matrix for the con¬ 
fining potentials of a harmonic or a square well potential 
assuming that the short-range spherically symmetric po¬ 
tential induces partial waves with angular and azimuthal 
quantum numbers (im) = [(00), (10), (1 ± 1)]. In addi¬ 
tion, connections between the current formalism and the 
LFT theory are discussed. Section IV is devoted in the 
discussion of several applications of the Schwinger varia¬ 
tional iF-matrix as well as its comparison with numerical 
simulations. Several appendices give some of the details 
of the derivations. Finally Section V summarizes and 
concludes our analysis. 

II. BASIC DEFINITIONS AND RELATIONS 

This section contains the conventions used in this pa¬ 
per for the different Hamiltonians with additive poten¬ 
tials, the scattering wave functions as well as the basic 
scattering parameters (e.g. the iF-matrix). In addition, 
the definitions of the corresponding Green’s functions are 
also included. Since there are many possible conventions, 
this will allow for unambiguous definitions of parameters. 

A. Hamiltonians 

The scope of this work is to treat the interaction of a 
particle with different potentials that have qualitatively 
different character, with the focus on the case where the 
particle’s motion is unbounded in one direction. The con¬ 
tinuum solutions can be obtained and used for the calcu¬ 
lation of scattering parameters, photoionization and/or 
photodetachment cross sections, etc. The full scattering 
information is encapsulated in the full Hamiltonian H 
which reads: 

H\^) = {Hf + V)\^)=E\^), (1) 

where Hf = j2^ is the free particle Hamiltonian and /i 

indicates its mass. V corresponds to the operator of the 
total potential, \tjj) is the scattering wave function and E 
is the total energy. 

In the following we assume that the operator V in the 
full Hamiltonian H is the sum of two potentials 14 and 
Vc, which possess simple functional forms in different co¬ 
ordinate systems. The two potentials are such that the 
wave function is easily calculated for each potential in¬ 
dividually. More specifically, the 14 will be a smooth 
potential that extends over large separation distances 
and affects the asymptotic boundary conditions in j-i/)). 
In this study, the V/ is considered to bound the mo¬ 
tion of the particle in all the degrees of freedom except 
one. For example, this potential could be a harmonic 
one, Vc[r) = [p/+ y'^) with frequency wj,, which 
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separates in cylindrical coordinates, or an infinite square 
well on the x — y plane which separates in Cartesian co¬ 
ordinates. The 14 (r) is a potential that is non-zero in 
a small volume r and could be a spherically symmetric 
potential at the origin. 

One key observation concerns the length scale separa¬ 
tion between the Vg and Vc potential terms. With this 
assumption in mind two new Hamiltonians are defined 
which use solely one potential term each. The He is the 
Hamiltonian that only contains Vf. and does not contain 
the W. 

Hem=E\^e), ( 2 ) 

where |'0c) is the corresponding wave function and E is 
the total energy. The solutions are analytic functions 
and are used in the following to define the asymptotic 
form of the scattering wave function. In addition, the 
corresponding collisional information can be expressed 
in terms of these states. 

Consider now the Hamiltonian Hg which contains only 
the Vg potential. Then the corresponding Schrddinger 
equation reads 


Hg 1^,) = E \i,g ), (3) 

where the eigenfunctions of this Hamiltonian have a sim¬ 
ple form because they do not include the complicated 
boundary conditions that arise from the confining poten¬ 
tial, Vc- Intuitively, Eq. m describes the motion of the 
particle at small distances where it experiences solely the 
short-range potential I 4 ■ 


B. Scattering wave functions and /f-matrix 

In this subsection the scattering wave functions of the 
Hg and He Hamiltonians are defined. More specifically, 
for the Hamiltonian Hg the corresponding potential term, 
namely I 4 , possesses spherical symmetry and is short- 
ranged. Therefore, the \tjjs) wave functions in Eq. ([3]) 
are expressed in spherical coordinates as radial functions 
times spherical harmonics. At distances r > tq, the Vg 
is negligible, namely 14(ro) ~ 0, therefore the AT-matrix 
normalized radial functions are phase-shifted from the 
free particle wave function according to the following re¬ 
lation: 

(r|'0s) = 

= \j-;^[u{kr) - tan(5^ ne.{kr)\Yim{^) (4) 

where the parameter k is defined from E = 
and the functions represent the spherical har¬ 

monic functions defined at angles H. The spherical Bessel 
[jiikr)] and Von Neumann [ni{kr)] functions correspond 
to the regular and irregular solutions of the Hamiltonian 
Hg. The term denotes the phase shift induced by the 
potential I 4 . 


In the following, the energy-normalized regular and ir¬ 
regular solutions of the He Hamiltonian are defined to¬ 
gether with the corresponding scattering matrices. The 
symbol tpe represents the regular functions at the origin 
and Xc represents the irregular functions. Since the I 4 
potential confines the motion in the x— and y—directions, 
then the parity in the z-direction, namely H^ is a good 
quantum number and hence the asymptotic form of the 
regular wave functions are expressed as even and odd 
solutions 


= +I 
= -l, 
(5) 

and the irregular wave functions possess the following 
form in terms of even and odd solutions 





sin(fcQ,|z|), for 
-|fyCOs(fca|2:|), 


Hz — -|-1 

for Hz = -1, 


( 6 ) 

where the 4>Q,(a;,j/) are a complete basis of orthonormal 
eigenfunctions of the transverse part of the He Hamilto¬ 
nian and their specific form is dictated by the particular 
type of confining potential 14 - The superscript notation 
(e, o) in Eqs. and (|ni) denotes the even (H^ = -1-1) and 
odd (H^ = —1) solutions, respectively. The total energy 
is given by the relation E = h^k^/2y + Ea where E^ 
is the energy of the transverse function These 

equations assume E > Ea assuring therefore that the so¬ 
lutions possess an oscillatory behavior. If the E < Ea, 
then there are functions that are real exponentials in z. 

The scattering solutions of the full Hamiltonian H at 
energy E can be expressed at small distances as linear 
combinations of the eigenfunctions of He Hamiltonian at 
this same energy E. Moreover, if the scattering potential 
is an even function of z, then the solutions separate into 
sums over the even functions or sums over the odd func¬ 
tions. In this situation, the asymptotic form of the exact 
wave function of the Hamiltonian H can be written as 




(7) 


where the elements K^^a'^ indicate the AT-matrix ele¬ 
ments of even or odd states in the z-direction. The a, j3 
indexes label the open channels {E > Ea,p) meaning that 
the corresponding wave function possess an oscillatory 
behavior as |z| —>■ 00 . On the other hand asymptoti¬ 
cally the closed channels {E < Ea) are described via the 
following functions 

T^®)(r) = <i>a{x,y) exp{-Ka\z\) ( 8 ) 

where the index a is here introduce in order to separate 
the closed from open channels. The energy of the closed 
channels is given by the relation E = — A^K^/2/i-|- Ea- 









4 


Having specified the scattering eigenfunction of the 
Hamiltonian, the corresponding i^-matrix fulfills the fol¬ 
lowing relation: 

Kap = -7r(^c,a|K|'0/3) (9) 

where a, /? label open channels. In addition, the 
I'i/'c), |'0) are the exact solutions of He and H Hamil¬ 
tonians, respectively. 

Equation ® can be expressed in a more general way 
by substituting the exact relation Vs = (E—Hc) IV'/?) 
and integrating by parts twice. Using this, the volume 
integral in Eq. (|9]) can be recast as a surface integral that 
is a simpler expression for the corresponding iiT-matrix: 

J Wz[^i%*{r),ilj^p\r)]dxdy, ( 10 ) 

where the term Wz[-] indicates the Wronskian with re¬ 
spect to the 2 -direction and is evaluated at large enough 
l^l that the closed functions of Eq. m are effectively 0. 
The closed channels are in the full wave function and 
could lead to resonances but do not contribute to the 
surface integrals at large 2 . 

C. Free particle and confining Green’s functions 


where p a vector lying into x — y plane. The index oq 
denotes the last open channel for a given total energy 
E. As with the Gj , the standing wave Gc can be ob¬ 
tained from the real part of G^. Note that this property 
does not hold in the case of a charged particle motion in 
magnetic fields. 

One important observation is that every Green’s func¬ 
tion diverges as ri 2 —>■ 0 in the same way for any energy 
not at a threshold ofHe'. 

Ganyiri,v,) = -^^^+0{r°,) ( 14 ) 

for small ri 2 ; this result also does not depend on whether 
the G obeys standing wave or incoming/out going bound¬ 
ary conditions. This means the difference of any two 
Green’s functions is finite everywhere except for energies 
exactly at threshold. In the expressions below, we will 
arrange to only have integrals that involve the differences 
of Green’s functions times functions that are finite and 
non-zero in a finite region. Thus, our expressions will au¬ 
tomatically be finite without the need for regularization 
procedures. 

In this work, the standing wave solutions are employed; 
therefore the corresponding Green’s functions, i.e. G/ 
and Gc denote the principal value ones. 


The Hamiltonians introduced in Eqs. m and dH) per¬ 
mit us to straightforwardly define the two corresponding 
Green’s functions. 


Gf^ 


1 

E-Hf 


Gc 


1 

E-He' 


( 11 ) 


where Gc only includes the confining potential, and G/ 
is the free particle Green’s function for the Hamiltonian 
with no potential energy. The specific asymptotic bound¬ 
ary conditions in Eq. (El) are determined by how the pole 
is handled. 

The outgoing/incoming free particle Green’s function 
G/ in the position representation can be written as 

~ II p±ife’'l 2 

(ri|G^|r 2 ) = G^(ri,r 2 ) =, ( 12 ) 


where ri 2 = |ri — r 2 | and the “-b” (“—”) denotes the out¬ 
going (incoming) Green’s function. The standing wave 
G/ is derived simply by taking the real part of either 

The outgoing/incoming Green’s function when only 
the confining potential is present reads 


(ri|G^|r 2 ) = G^(ri,r 2 ) 

cuo 


= E 


Q!=l 

00 

E 

ao + 1 


h'^kr. 


$^(Pl)$a(P 2 )e 


:Lika\zi — Z 2 \ 


h?K, 


-$a(pi)$*(p 2 )e 


— Koi\zi—Z2 
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III. SCHWINGER VARIATIONAL PRINCIPLE 

This section presents a derivation of the AT-matrix in 
terms of the Schwinger variational principle. The main 
focus of this derivation is the standing wave solution and 
the corresponding Green’s function. However, note that 
a similar derivation can be carried out in cases of incom¬ 
ing/outgoing boundary conditions. 

Using the Lippmann-Schwinger equation, \'4>a) = 

|'0c,a)+GcUs I'i/’a), and its hermitian conjugate in Eq. (O, 
the following identity for the AT-matrices is fulfilled 

Kap = -TT {'lpc,a\Vs\i’p) = -TT | K | V’c./s) , (15) 

where the relation G\ = Gc is used. 

A Schwinger-type variational expression for the K- 
matrix can be obtained with the help of Eq. ED and 
the trial functions = \ip) -b \6ip) of the full Hamil¬ 

tonian H. 

+ (^W|IG GcV;|<^)], (16) 

where the variational expression for the A"-matrix, K™, 
equals the exact A"-matrix, i.e. Eq. plus terms of 
order 5ijP". 

EZ" = - WalK - VsGeVs\5i^p)] 

= Ke,p + OiSijZ (17) 
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The trial functions are often written as linear com¬ 

binations of a basis set of functions \yj). 

3 

Substituting Eq. (IT51) in Eq. (fTBl) . the coefficients Cja 
are specified by requiring Eq. (HU) to be variationally 
stable, i.e. dK'^^/dCjp = 0. This yields the basis ex¬ 
pansion version of the Schwinger variational expression, 
which reads 

^^cAVs\y,) {yy\VsH,,p) (19) 

33 ' 

where Mjij = {yj'\Vs — VsGcVslyj). In addition, the trial 

functions \'tpa'^) can be evaluated via the coefficients Caj 
which fulfill the following relation 

Cjp = {y3'\t\M ■ ( 20 ) 

f 

One of the main points of a variational principle is that 
an exact i^-matrix is obtained if the |'!/'*^‘^)’s equal the 
exact {ipYs. For the Schwinger variational principle, all 
the corresponding expressions for the if-matrix involve 
terms of 14 This implies that the exact result for 

the if-matrix is obtained for the Schwinger variational 
principle even if the satisfies the somewhat looser 

condition 14(r)(5^(r) = 0 which will be exploited in the 
next section. 

In particular this condition can be satisfied at large 
and small distances r. At large distances, the potential 
14 vanishes implying that the 5'0(r) parts of the trial 
functions |^^*^) which are nonzero in this region do not 
contribute to the if-matrix calculation. On the other 
hand, the trial functions at small distances can be chosen 
such that they are exact solutions of the full Hamiltonian 
H yielding therefore (5'0(r) = 0 at small r. This imposes 
a constraint on the choice of the trial functions or the 
basis expansions that is employed in the next section. 
With these choices, the exact AT-matrix is obtained even 
if the is not an eigenfunction of the full Hamiltonian 
H over all space. 

An additional feature is that if 14 (r) 7 ^ 0 in a small 
region of space (and for a limited range of low energies), 
only the low-£ angular momentum partial waves of 
are needed to calculate the if-matrix. This holds only 
away from high-£ resonances due to the fact that as the 
angular momentum (. increases the amplitude of the cor¬ 
responding wave functions in the region of 14 (r) 7 ^ 0 van¬ 
ishes. Therefore, these high-£ states will not contribute 
to the scattering information. 


IV. K-MATRIX FOR TWO POTENTIALS: 

APPLICATION OF THE SCHWINGER 
VARIATIONAL PRINCIPLE 

An important special case is when the potential l4 = 
H — He is nonzero over a range small enough that the 
potentials in He Hamiltonian are considered effectively 
zero (or constant) over the range of I 4 . This will allow 
us to simplify all of the integrals in Eqs. (HU) and (HOI) 
and in some cases the corresponding results are simple 
analytic functions. To reduce the amount of new mate¬ 
rial, we have restricted the scattering systems to have no 
z-dependence to the potential and for the confining po¬ 
tential to have symmetry in x, y. Finally, an important 
consideration for this section is to obtain expressions for 
the M-matrix elements, Eq. dni), as a simple analytic 
term and a term involving an integral over finite func¬ 
tions. 


A. General expression of the K-matrix for 
short-ranged spherical symmetric potential 14 

As discussed at the end of Sec. uni only low-^ angu¬ 
lar momentum partial waves contribute to the calcula¬ 
tion of the AT-matrix when Vs is assumed to be short 
ranged. Therefore, for such a short range potential only 
the first two angular momentum states will be consid¬ 
ered, i.e. £ = 0 and 1. This approximation will allow 
us to obtain closed form expressions of all the A'-matrix 
elements in terms of the £ = 0 and 1 phase shifts. 

Since the full potential V is symmetric in z, the K- 
matrix separates into even- and odd-channel blocks. For 
each parity state of the He Hamiltonian, only one basis 
function is used in Eq. m which is the solution of H 
Hamiltonian. This will provide the exact scattering infor¬ 
mation given that the length scale of the potential terms 
14 and Vc in the full Hamiltonian H are well separated. 
In other words, the length scale separation implies that 
in the region of I 4 ^ 0 the states |'0s) (see Eq. ([3])) will 
not be modified by the confining potential 14- 

To be precise, if there is a potential that only scat¬ 
ters s- and p-waves at some specified (low) energy, then 
the potential only extends over a very short range. The 
wave function in the region near the scatterer, at posi¬ 
tion Tg, is most easily represented in spherical coordi¬ 
nates. The difficulty will be to evaluate the integrals 
involving the \tjjc) and the corresponding Gc which will 
often be most easily represented in a different coordinate 
system. The reason that the exact scattering informa¬ 
tion is obtained from the Schwinger variational principle 
is that the exact wave function of the full Hamiltonian 
(including large r boundary conditions) must have the 
form '0(r) = C'0s(r) + C('^)- Note that Y’s is the spher¬ 
ical wave centered at Vg that does not obey the large-r 
boundary conditions of the full Hamiltonian. The func¬ 
tion ^ = '0 — Cil^s only has £>2 partial waves at which 
implies I4(r)^(r) = 0 everywhere. This means the exact 
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A'-matrix can be written as 


|'*/’s) ('0s IV'c./s) 

±\.(yQ - TT ,v ^ .V ^ 5 (^i) 

{i^s\Vs-VsG,Vs\^s) 

where each term is now a single matrix element. The 
terms in the numerator give a finite result because the 
Vs I'0s) is non-zero only over a finite region. Although the 
term in the denominator is not obviously finite because 
the Gc diverges in proportion to l/|ri—r 2 |, the derivation 
below shows that this term can be written as convergent 
sums. 

The evaluation of the term in the denominator of 
Eq. (| 2 T]) is complicated by the fact that the Gc is simple 
in a different coordinate system from that of |0s) and 
14 [H, [2^-[33 ■ However, these complications can be 

avoided by adding and subtracting the Green’s function 
for a free particle, G/, 

(0014 - 14Gcl4|0s) =Is + Ds (22) 
with the new parameters A and Dg defined as 


nonzero only over a small region, the integral should only 
be performed over a volume r that contains this region. 

The Appendix |A] gives the prescription for evaluating 
these matrix elements for £ < 1. The basic idea is to 
convert the volume integral into a surface integral which 
is possible due to the length scale separation between the 
range of 14 and when the 14 becomes non-negligible. The 
main results are: 


(0c!al^^*l0^<.oo) = (^-^tanSo^ (26) 




^(-Itan.,) (27) 




127r / 1 


where (^±$ 0 ) = {d^a{^)/dx^id^a{j^)/dy)/2. 


(28) 


Is = (0014 - 14G/14I0,) = --tanSi (23) 

TT 


C. Analytic expressions for the Ds matrix elements 


where the | 0 s) indicates a state of a particular angular 
momentum £ possessing the form of Eq. o and 

Ds = (0014AG1/000 = (0014(G/ - G,)14|00 (24) 

These manipulations have the advantage that the Is is 
simply related to the no confinement phase shift and the 
G/(ri, r 2 ) —Gc(ri, r 2 ) function is finite for all values of ri 
and r 2 . Thus, the expression in Eq. (1^ is finite because 
the 14 10s) is finite with a finite range and its square is 
multiplying a finite function. Note that a similar differ¬ 
ence of Green’s functions appears also in the methods 
presented in Refs. iHl yielding thus finite results. 

One important observation is that the AT-matrix in 
Eq. (I?T]) is essentially separated in factors of matrix el¬ 
ements. This is crucial since it allows us to treat each 
factor separately in a pedagogical manner. The follow¬ 
ing subsections are devoted to providing their explicit 
analytical expressions. 


B. Matrix elements with 14 |0s) 

It is evident from Eq. (l2T]l that most of the integrals 
involve terms of 14 10 s); thus this subsection mainly focus 
on the integrals of the form 

(0c|H00s) = J 0c(r)14(r)0s(r)d^r, (25) 

where r is the smallest volume that contains the region 
where 14(r) 7 ^ 0. Recall that the function 0c(r) is an 
eigensolution of the He Hamiltonian. Since the 14 (r) is 


This subsection focuses on the Dg matrix elements. 
One of the main assumptions here is that the 14 potential 
is centered symmetrically with respect to the confining 
potential; for example, 14 is at the origin and Vc{x,y) = 
Vc{y,x) and Vc{x,y) = Vc{—x,y) etc. This allows us to 
resolve possible challenges with respect to the scattering 
of £ = 1 , \m\ = 1 . 

The Dg matrix elements for £ = 0, m = 0 and for 
£ = 1 , m can be written in the compact form 


Dg^im — 



tan 5(, 
[a_Lfc]2^+i 


mm{E) 


(29) 


where the a_L is a relevant length scale for the con¬ 
finement potential I 4 and ? 7 fm(A)-function is defined in 
terms of the AG defined in Appendix[Bj For £ = 0, m = 0 
and £ = 1, 'm, they are defined as 

yooiE) = ^aj.AG(0) 

7?io(A) = -^aiAG,,(0) 

m±i{E) = -^aiAGpp(O) (30) 

It is important to note that the expression for AG 
converges for non-zero Az without the need for regular¬ 
ization. In order to investigate in detail the Az depen¬ 
dence at Ap = -^ <C I, the function r(ri,r 2 ,A) = 
27rfi^a_LAG(ri, r 2 )/// is defined which reads 


r(ri,r 2 . A) 


27r^ai$a(pi)4>;(p2) 

CX. — 1 


sin(A:Q,|Az|) 

ka 
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2^ ai®a(Pl)®a(P2) 

a=Q:o + l 

C0s(fcfi2) 

ri 2 


3 —Kq I Az\ 


(31) 


where the Ea^ < E < Ea^+i defines the open and closed 
channels. Parameters with a tilde have been scaled by the 
length scale a_L to become dimensionless: Az = Az/a±, 
k = ka±, ka = kaa±, and ka = Kaa±. In the expres¬ 
sion for r, the first term is the contribution from the 
open channels to Gc, the second term is the contribution 
from the closed channels to Gc, and the last term is the 
contribution from Gf. 

The three rj^miE) (see Eq. (1501) ') functions are derived 
from the r(ri, r 2 , E) in the limit Az —>■ O’*" and Ap 1. 
Physical considerations from Eq. m means the expan¬ 
sion of the r(fi, f 2 , E) will have the form 


r(ri,r2,E) ~ pooiE) + ziZ2 


+ 


miiE) 


3 

m-i{E) 


-Pl ■ P2- 


(32) 


The physical meaning of the pimiE) terms is that they 
encapsulate the impact of the confining geometry, i.e. 14, 
on the short-range spherical symmetric potential 14 ■ In¬ 
deed, Pirn conveys the information of each spherical wave 
of £m character induced by the 14 at short distances is 
distributed over the asymptotic a-channels imposed by 
14 at long distances. Evidently, the explicit form of the 
Pirn terms depends on the specific type of the confining 
geometry. Therefore, in the Appendix [Cl the pim{E) 
terms are evaluated in detail for two types of confining 
potentials, e.g. the harmonic and the square well poten¬ 
tial. 


D. Test of Taylor series expansion of AG 



Figure 1. (color online) The quantity F of Eq. (1321) as a func¬ 
tion of Az/a±ior the confining potential of the harmonic os¬ 
cillator infinitesimally above the threshold energy E = /kJx 
where pi = p 2 = 0. The black solid line is denotes the sum 
in Eq. (IC3t . The red dashed line depicts the approximation 
Eq. (1321) using Eq. (HU). The blue dotted line indicates a nu¬ 
merical extrapolation of r(ri,r 2 ,i?) = r(Ai —>■ 0’^,Ap —>■ 
0, E) ~ [4F(A5,0, E) - r(2A5, 0, E)]/3. 



Figure |T] explores the various approximations on the 
quantity T which were discussed in Eqs. dSD and ()C3I) 
considering that the confining potential is a harmonic 
one. The solid line in Fig. |T] depicts T—quantity from 
Eq. (ISTI) at energy infinitesimally above the threshold 
E = huj± as a function of Az for pi = p 2 = 0. The 
latter means that the azimuthal quantum number m is 
set equal to zero. The red dashed line in Fig. |T] cor¬ 
responds to Eq. (IC3I) which is obtained by calculating 
the poq{E) and piq{E) using the method described in 
the appendix |C] The red dashed and the black solid line 
are in excellent agreement, especially for small Az. This 
highlights the validity of the approximations that were 
considered in Eq. (IC3I) . In addition, the blue dotted line 
refers to pi=o,m=o{E) parameter which is defined as a 
linear combination of the r(ri,r 2 ,E) functions cancel¬ 
ing in this manner terms of the order |Azp according to 
the relation r(ri,r 2 ,£') = r(Az — >• 0+,Ap —>• 0,E) ~ 
[4r(A5,0, E) - r{2A2, 0, E)]/3. Therefore, in Fig. H] the 
dotted, solid and dashed lines converge to the same finite 


Figure 2. (color online) The quantity F of Eq. (1321) for the 
square well infinitesimally above the threshold energy E = 
ft^7r^/[^ax] as function of AS for pi — p 2 = 0. The black 
solid line denotes the sum in Eq. (IC7l) . The red dashed line 
indicates the approximation Eq. (EH using Eq. dig. The 
blue dotted line is a numerical extrapolation of r(ri, r 2 , E) = 
r(AS ^ 0+, Ap 0, A) ~ [4r(AS, 0, E) - r(2AS, 0, E)]/3. 


value of pi=o,m=o without exhibiting any divergences. A 
similar analysis can be carried out for the pi=i^rn=i{E) 
and pe=i,m=±i{E), and it shows excellent convergence 
between Eqs. m and (IC3I) . 

Similar to Fig. |T1 Fig. |2| investigates the validity of 
Eqs. (1511) and (IC7I) for the case of a square well con¬ 
fining potential. Fig. [5] mainly depicts the quantity T 
as function of |AS| for pi = P 2 = 0 infinitesimally 
above the threshold energy E = h?Tr‘^/[pa^]. The black 
solid line indicates Eq. (EH) and the dashed line corre¬ 
sponds to Eq. (IC7I) . It is apparent that both relations 










are in excellent agreement in particular for Az —>■ 0 
suggesting in this manner the regime of validity of the 
Eq. (IC7I1 . In addition, note that the black solid and the 
red dashed line converge to the same value as Az —t 0, 
i.e. 7700 (E). Furthermore, the blue dotted line in Fig. 
consists of a linear combination of the r(ri, r 2 , E) quan¬ 
tity such that terms of order |Azp defined by the re¬ 
lation F(ri,r 2 ,E) = F(Az —^ 0+,Ap —>• 0,E) ~ 
[4r(Az, 0, E) — r{2Az, 0, E)]/3. Evidently, in Figl2] this 
linear combination (blue dotted line) is constant in the 
interval of small Az and corresponds in essence to the 
value of 7700 (E). 


V. SCATTERING OBSERVABLES IN TERMS 
OF THE SCHWINGER VARIATIONAL 
A-MATRIX 


This section utilizes key formulas developed in Section 
nYi in order to explicitly connect the Schwinger varia¬ 
tional ET-matrix of a two potential Hamiltonian with all 
the relevant scattering observables such as the photoab¬ 
sorption cross-section. In addition, the transmission and 
reflection coefficients for quasi-one dimensional Hamilto¬ 
nians are derived in terms of the Schwinger variational 
ET-matrix. This derivations permits us also to discuss 
the connection of the present study with the theoretical 
framework of the local frame transformation (LFT) the¬ 
ory. We reiterate that the derivation is for the case of the 
confining potential being independent of z and symmet¬ 
ric in xy, examples that do not follow these restrictions 
will be presented elsewhere. 


A. The explicit analytical expressions of the 

A-matrix and its connection with the LFT theory 

In this subsection the formulas derived in the subsec¬ 
tions IIV HI and II V Cl are the main constituents of the 
Schwinger variational A-matrix in Eq. (EH- Recall that 
the short-range scatterer is assumed to be placed at a 
symmetry point of the confining potential. This implies 
that the A-matrix will be non-zero in blocks depending 
on whether the asymptotic channel a couple to £m states 
or not. Note that each a-channel is described either by 
even or odd function in z-direction, thus these states cou¬ 
ple to im spherically symmetric states if -I- tti is either 
even or odd, respectively. 

In this manner, substituting Eqs. (l23l) . (l26l) and (l29l) in 
Eq. (I2ID the A-matrix is provided for a given im pair of 
quantum numbers, where the type of confining geometry 
is not specified. The latter permits us to provide a generic 
form of scattering matrix elements which read 


7 y 


;2^+l 


(E) 


1 + af+\E)r,e^{E) 


[Wt] 




(33) 


where the superscripts (e, o) refer to even or odd states in 
z-direction, respectively. Note that the terms a^^“''^(E) = 
a^^+\E)/af+^ where af+\E) = - tan6e{E)/k^‘^+^ 
correspond to i-wave energy-dependent scattering length 
(i = Ojor volume (i = 1). It is important to note that 
for some long range fields, ie. the polarization poten¬ 
tial, the scattering volume diverges in the limit of zero 
total energy. On the other hand, in the problems of the 
present study the additional confining potential sets a 
lower limit in the total energy other than zero avoiding 
in this manner the divergence of the scattering volume. 
However, in the regime of extremely weak confinement 
the lower bound in the total energy tends to zero yield¬ 
ing eventually a divergent scattering volume which must 
be analyzed more carefully. r)im{E) is the next energy de¬ 
pendent parameter. For specific im state pirniE) param¬ 
eter intrinsically possesses threshold singularities which 
occur below (a)-channel thresholds (see FigsE] and H]). 
The matrix elements Ua,im or its conjugate transpose 
capture the coupling of the short- and long-range 
physics. In other words, the matrix U contains the infor¬ 
mation that relates the short-range quantum numbers i 
and m with the a ones which are fully dictated from the 
confining potential, i.e.. long range physics. 




27ra 


^<l>„(0,0) 

T 


y'67ra^fca$a(0,0) 

«•-)!( 


£ = m = 0 
£ = 1 , m = 0 
£ = 1, m = ±I 


(34) 


where the operators are defined below Eq. (I28p . 

At this point it becomes evident that Eq. (IMl) is in 
essence the local frame transformation (LFT) which per¬ 
mits the the projection of the a—states on the £, m spher¬ 
ically symmetric states. Note that U in Eq. (IM)) and the 
LFT transformation U in Refs. @,0 are connected ac¬ 
cording to the relation = Ua,im{ka±Y~^^^^. There¬ 
fore, specifying the type of the confining potential, i.e. 
harmonic oscillator or square well, one obtains the local 
frame transformation either from cylindrical to spherical 
coordinates Q or from Cartesian to spherical coordinates 
[T^ . In particular, in Ref. [l^, the i = m = 0 derived 
LFT is identical to Eq. (IMl) . Similarly, the LFT of Ref. @ 
provides the same results as Eq. OMl). However, for com¬ 
pleteness reasons it should be noted that there is a mis¬ 
take in Eq. (27) of Ref. Q, because strictly following the 
notation of Ref. @ would yield the new formula, namely 

Uni{m) = {7,^y'’(77i)(2a)^/^[7r!/(7r-|-|77^|)!]“^/^[i(27^,-|-r?^-|- 

|77^| -I- 

In addition, the A-matrix equation, namely Eq. dSSl), 
is exactly the same as the corresponding physical K- 
matrices of Refs. [13,Hill- However, there is a concep¬ 
tual difference between the framework presented here and 
the framework in Refs. PHIS- This difference arises 
from the fact that the A-matrices in Refs. PH [3 did 
not obey the physical boundary conditions for the closed 
channels; therefore a closed channel elimination was em¬ 
ployed in order to obtain the physical A-matrix which 
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obeys the proper boundary conditions. This procedure 
yield ri£^[E) functions which involve only the sum parts 
of Eqs. (IC3I1 and (IC7I1 . Thus, the TFT leads to diver¬ 
gences. This singular behavior was removed in a sec¬ 
ondary step that required thoughtful implementation of 
auxiliary techniques of regularization , such as Riemann 
zeta function regularization [T^ [l^ or residue regular¬ 
ization [l^ . On the other hand in the present theoret¬ 
ical framework no such techniques need to be employed 
since by definition rji^{E) functions are finite. However, 
note that r]im{E) possess intrinsic threshold singularities 
at energies infinitesimally below the channel thresholds 
of the confining potential. Furthermore, the Schwinger 
variational formalism provides an intuitive understand¬ 
ing of the regularization techniques. As is shown in the 
Appendix [Cl the ijim (E) functions consist of a difference 
of a sum and its integral representation, whereas the lat¬ 
ter originates from the free space Green’s function. This 
piece of information is absent from the physical AT-matrix 
in the TFT theory and it is incorporated by various reg¬ 
ularization schemes. Perhaps a comparison between the 
derivation of the LFT and the present results could lead 
to a deeper understanding of the source of the divergences 
in the LFT, and deserves to be studied in the future. 


B. Photoabsorption Cross-section 


In order to obtain photoabsorption cross sections, we 
define the dipole matrix elements which describe the 
transitions from the initial state to the incoming wave fi¬ 
nal state that possesses only outgoing waves in the a — th 
channel yielding the following expression: 

(37) 

/3 


where the iC-matrix has rank one for each im pair, with 
the relevant eigenvector being proportional to the Da- 
This allows a relatively simple expression for the dipole 
matrix elements. 




Dt.m{E) 


l + af+\E){l^UE)-^N^tmy 


(38) 


with 


iVL = E (39) 

a—0 

where the summation index a indicates a sum over only 
open channels. Then the total probability is proportional 
to \D~\^ = \^a \ which is the sum over the partial 

absorption terms Pa = The sum reduces to the 

relation 


Since the Schwinger variational A'-matrix is defined 
one can also define the electric dipole matrix elements 
and consequently the photoabsorption cross section. For 
example, these relations are needed for the treatment of 
negative ion photodetachment in the present of external 
fields. Initially, the expression for the trial function coef¬ 
ficients on the basis of the \tjjs), he.. Fq. (EUl) are defined 
as 


Ci 


m,a. 


l + af+\E)r,e^{E)' 


(35) 


Since we are interested in computing the cross section 
that describes, for example, the excitation of a negative 
ion by photon absorption, the dipole matrix elements at 
small distances are defined by the relation Dim{E) = 
ii’Cmls ■ Note that the term e • f is the dipole 

operator with i denoting the polarization vector and the 
state ipinit corresponds to the initial state of the negative 
ion. Having defined the dipole matrix elements at short 
distances, one can derive the dipole matrix elements that 
describe the transition from the initial state to the a — th 
state which includes the confining potential. Therefore, 
the dipole matrix elements read 


Ea — E^^i^E^C^^a — 


DUE)Ul^,a 

I + af+\E)ini^{E) 


(36) 


where the element Dim{E) = Egm{E)/{kaj_Y~^^E typi¬ 
cally has less energy dependence than Dim{E). 


\E-\'^ = 


NLDUE) 


[1 + ~af+\E)rjUEW + [«r^(7^)^^L]^ 


(40) 


C. Transmission and reflection coefficients 

This subsection contains expressions for the transmis¬ 
sion and reflection coefficients of quasi-one dimensional 
Hamiltonians in terms of the Schwinger variational AT- 
matrix in Eq. dSD- This situation models the scattering 
aspects of two-body collisions in the presence of a trans¬ 
verse harmonic potential or else for collisional complexes 
of light particles being scattered by heavy ones in the 
presence of a confining transverse potential. 

By superposing the even and odd wave functions, the 
solution with an incoming wave in the a — th channel can 
be obtained when s- and/or p-wave scattering dominates. 

+ ( 41 ) 

P 

where the (r) only has a wave incoming from negative 
2 in channel a. In the case of identical particle collisions 
having only s-wave (p-wave) character, the odd (even) 
AT-matrix will vanish, namely = 0 (AT^®) = 0). On 
the other hand in the case of distinguishable particles, s- 
and p-waves could potentially both be present, and then 
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Figure 3. The 7700 (black solid line), rjio (blue dotted line) and 
7711/20 (red dashed line) for harmonic oscillator confinement 
as a function of the scaled energy e = E/{hLo±). 


both even and odd parity if-matrices can contribute in 
the 7 / 1 + (r) wave function. 

From the preceding equation one can define the ma¬ 
trices of the corresponding reflection and transmission 
amplitudes which fulfill the relations 


iK ^ 1 

1-iK * ~ l-iK' 


(42) 


In the case of s-wave scattering the reflection and 
transmission amplitudes can be expressed in terms of the 
energy dependent scattering length and the matrix lA : 

_ 7 , iao{E) . 

— ~^a,00 1 , - /rnw - lT^\ i ■ Ar2 1 ^ J00,/3 

l + ao(^)(?7oo(-E^)+^A^oo) 

~ ^I3a H" '^13a 

Note that similar expressions can be obtained for p- and 
s -f p-wave scattering but it is straightforward to derive 
them if needed. 

From these matrix elements, the scattering probabil¬ 
ities can be obtained. The probability for an incoming 
wave in channel a to reflect into channel (3 is p while 
the probability to transmit into channel [3 is Itapp. Note 
that the reflection and transmission probabilities are the 
same except when (3 = a. This is understandable on 
physical grounds: a zero range potential scatters equally 
into +z and — z but the f3 = a channel has different 
interference with the incoming wave in the +z and — z 
directions. 


VI. RESULTS AND DISCUSSION 

This section focuses on the various aspects and insights 
provided by the present framework of the Schwinger vari¬ 
ational /f-matrix approach and its comparison with the 
LFT results. In addition, a set of selected applications of 
the current approach is thoroughly discussed to highlight 
its robustness and clarity. A first test bed is taken to be 


the study of photodetachment processes in the presence 
of an external confining field where the closed channel 
physics is taken into account, thus going beyond previ¬ 
ous studies 

Recall that the closed channel physics is fully encapsu¬ 
lated in the 77 ^^ (A) functions which depend on the par¬ 
ticular im states of the short range potential and on the 
specific type of the confining potential. 

In other words, r]irniE) corresponds to a collective pa¬ 
rameter, which emerges from the strong mutual coupling 
of all the closed channels caused by the spherically sym¬ 
metric short range potential at small distances. In terms 
of the LFT theory, this collective parameter is described 
by the eigenvalues of the closed-closed partition of the K- 
matrix and, as was mentioned above, these AT-matrices 
are of rank 1 despite the fact that their dimension is 
M X M with M —>■ 00 . This particular property ensures 
that there exists only one nonzero eigenvalue representing 
the effect of all of the closed channels. Incorporating this 
concept in photodetachment processes yields non-trivial 
resonant features; therefore, in the following subsection 
the energy dependence and broad features of the rjemiE) 
functions are studied for various confining potentials. 


A. Behavior of rjirniE) 


1. Harmonic and infinite square well confining potentials 


Figures [3] and m depict the quantity riem{E) as a func¬ 
tion of the dimensionless energy e for the harmonic 
and infinite square well confining potentials, respectively. 
Note that the scaled energies e fulfill the relations e = 
E/huj± and e = E/lAh^Tr"^/[2p,a\]) for the harmonic (Ap¬ 
pendix 1C II) and square well (Appendix 1C 21) potentials, 
respectively. 

The expressions for r]im for the harmonic and square 
well potential are given in Eqs. (IC4I) and (IC8I) . respec¬ 
tively. These parameters are expressed in terms of sums 
and integrals involving the scaled momenta and Az which 
can be simplified to show only the dependences on the 
quantum numbers. For the harmonic oscillator, Eq. (IC4I) 
becomes 


mo = 


g—2-\/n— l'oA5 /*oo ^—2-\/n—i^oA.z 


n>L'o 


\/n - r'o 


y/n - Vo 


dn 


^ — 2y/n—i'oAz 


n>vo 

nOO 


7710 = -12 ^ a/uTUi/o e 

(n -k i)e- 2 V"-i'iAi 


-|-12 / ^/n — vq e 


m±i 


= E 


n>y-f 


\/n - ixi 
(n-k 

y/n - 1/1 


dn 


(44) 
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where the limit Az —> O'*' is understood in these expres¬ 
sions. nUm is given hy = (e — 1 — \m\)/2 and it relates 
with momenta in Eq. (IC4|) according to the rela¬ 
tion kn,m = 2i/n — Vm- For the square well confinement, 
Eq. (IC 8 |) becomes 


g- 27 r.y/ag-eAz eoo g- 27 r-eAi ^ 

= -rT=-4/ - - -ada 

Vag-e Va^-e ^ 


»7io = 






pCO 

-|-487r^ / \JcP- — 

JV^ 


' —ada 


m±i 


= 48-^E 


{a^ -f i)2e-2^\/^Ai 

\/a?-e 


.02/’“ TT 3 , 

— 487 r^ / - , - -a^da 


(45) 


where where the limit Az —> O’*' is understood in these 
expressions, is given by = (a^, -I- [|to| -I- l]/2)^ -I- 
(tty -I- 1/2)^ with Qfx, ay independently being 0,1,2,.... 

is related with the momenta kam in Eq. (IC8I) accord¬ 
ing to the expression Ka,m = 27rY^ a^ — e. In the inte¬ 
grals, the double integrals daxduy have been converted 
to polar coordinates and the integral over angle from 0 
to 7r/2 has been carried out. 

Figure [3] shows the results of Eq. (I44p for harmonic 
oscillator confinement for {£m) = [( 00 ), ( 10 ), (11)] states. 

We identify that the parameters 7700 (E) and 7710 (E) 
from Eq. (1441) correspond to the Hurwitz Riemann func¬ 
tions C(l/2,3/2 - E/2hu}±) and ({-1/2, 3/2 - E/2hu;±), 
respectively. Note that exactly the same formulas was 
shown in Refs. [13, [il. Furthermore, as was shown also 
in Ref. [s^, the 7700 (E) (see black solid line in Figl3]) 
and 7710 (E) (see blue dotted line in Figl3]) are periodic 
functions of e = E/{huj): rjirniy -|- 2) = 77^m(e)- 

However, the 7711 (E) function (see red dashed line in 
FigU) starts at a higher threshold energy because the 
lowest transversal state for tti = 1 is an hu) higher in en¬ 
ergy than for m = 0 and it is not periodic in e. Actually, 
the 7711 (E) consists of a sum of the periodic functions 
7700 (E) and 7710 (E). However, the corresponding 7700 (E) 
term is proportional to the energy, which explains the 
non-periodic character of 7711 (E). 

The 7700 and 771^1 functions diverge as the energy ap¬ 
proaches a channel from below due to threshold singulari¬ 
ties whereas the 7710 remains finite at each channel thresh¬ 
old. More specifically, the 7700 and 771^1 diverge due to the 
terms of the 1 /kq in the summation whereas the 7710 does 
not exhibit such features at threshold because the sum¬ 
mation involves terms proportional to Kq. Furthermore, 
at each threshold, the value 7710 = 2.494632... which is 
the same as in Refs. [l3, [T^. Similarly, the 7700 parame¬ 
ter acquires the value 7700 = —1.460354... by approaching 
the thresholds from above which is the same as in Refs. 
Em. The 7711 function has the values 7711 = —10.00944 



Figure 4 . The 7700 (black solid line), 7710/5O (blue dot¬ 
ted line) and 7711 /300 (red dashed line) for infinite square 
well confinement as a function of the scaled energy e = 
E/{Ah‘-K‘^/\2iia\]). 


at the lowest threshold and —18.77157 at the first ex¬ 
cited threshold when approaching those thresholds from 
above. On the other hand, when the threshold is ap¬ 
proached from below the 7711 function diverges and this 
behavior arises from the singular behavior of the 7700 . 

Figure m illustrates the r]im{E) functions for the infi¬ 
nite square well potential where [Im) = [(00), (10), (11)] 
states are considered. More specifically, the functions 7700 
(black solid line), 7710 / 5 O (blue dotted line) and 7711 /300 
(red dashed line) are shown versus the scaled energy 
e = E/( 4 fi^ 7 r^/[ 2 / 7 a^]). In contrast to Fig. |3l none of the 
77 ’s are periodic functions of e. As with Fig.O the 7700 and 
7711 go to 00 as the thresholds are approached from below 
whereas they acquire a finite negative value as the thresh¬ 
olds are approached from above. This behavior arises 
from a similar analysis to that for Fig. |3] but now involv¬ 
ing the individual terms in Eq. (1451) . As in the harmonic 
oscillator example, 7711 starts at a higher threshold energy 
because the transverse functions must have a node to con¬ 
tribute, since the azimuthal quantum number is m = 1 . 
Our threshold value for 7700 = —5.853459... is in excellent 
agreement with Ref. [l^ result from the R-matrix calcu¬ 
lation but is ^ 0 . 2 % different from the result of their reg¬ 
ularized analytical summation. In addition, from Fig. |4] 
the threshold values of the functions 7710 and 7711 acquire 
the values 7710 = 132.098... and 7711 = —642.613..., re¬ 
spectively. Note that 7710 has been divided by a factor of 
50 and 7711 by a factor of 300 to plot them on the same 
graph. This indicates that the effect of the confinement 
on p-wave scattering will typically be larger for the infi¬ 
nite square well than for the harmonic oscillator. 


2. Infinite Rectangular Well Confinement 

Due to the simplicity of the current formalism one can 
explore the properties of various confining geometries. In 
this section, the case of a rectangular well confinement is 
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Figure 5. (color online) The threshold value of r;oo (black solid 
line) and 7710/10 (blue dotted line) for the infinite rectangular 
well as a function of the ratio of conhning length scales, c 
dehned as |a;| < a_L/2 and \y\ < ax/(2c). Note 7700 goes from 
negative to positive at c ~ 4.89. 


considered. This particular geometry permits us to study 
the impact of lifting the degeneracy of the square-well 
potential. Furthermore, in the following the discussion 
is restricted on s-wave scattering which couples only on 
even z-parity states of the rectangular well and/or on p- 
wave collisions which couples odd z-parity states of the 
rectangular well. 

The rectangular potential well is defined as follows: 
V(x,y) = 0 when |a;| < ax/2 and \y\ < a±/(2c) and 
V(x, y) = 00 otherwise. Note that the parameter c con¬ 
trols the degree of asymmetry in the potential. The 
transverse eigenfunctions of the corresponding confining 
Hamiltonian He are well known and close to the origin 
obey the relation 

Of — 

^ V (46) 

a±^/c 

for all {axjCty) that are non-zero at the origin. The cor¬ 
responding eigenenergies for the states that are nonzero 
at the origin can be written as 

+ K/c]^) (47) 


where Ux = Soj, -|- 1 and Uy = 2ay + 1 where the cty 
can independently be 0,1,2, 3,.... 

Since we are interested in s— and p-wave scattering 
which respectively couple to even and odd-parity states 
of the rectangular well the prescription given in the ap¬ 
pendix 1C 21 is followed in order to obtain the correspond¬ 
ing parameters. Namely, the 7700 and 7710 parameters ful¬ 
fill the following relations: 


Voo 


4 

c 


E 

Qg>e 


g-27r.y/ag-eAz 


- 4 



— eAz 

\J o? — t 


—ada 


Figure 6. (color online) The 7700 parameter is depicted in pan¬ 
els (a) and (b) for c = 5 and c = 40, respectively. The 7710 /IO 
parameter is depicted in panel (c) for c = 5 (black solid line), 
and c = 40 (red dashed line). All the eto-parameters are plot¬ 
ted as a function of the scaled energy e = iJ/(4fi^7r^/[2/raxl)- 
The magenta arrows in panel (b) show the thresholds for the 
tight confinement in one direction, ie c 1, which acquire 
the values e = 1/4, 9/4, 25/4, ... 


mo 


AO 2 ^ _ 

ag>e 

fOO _ /TT 

+487r2 / 


(48) 


where where the limit AS —>■ 0+ is understood in these 
expressions and Oq obeys the relation Oq = (ax + 1/2)^ -I- 
(ay-l-1/2)^. In the integrals, the double integrals daxday 
have been converted to polar coordinates and the inte¬ 
gral over angle from 0 to 7 r /2 has been carried out. Our 
numerical tests demonstrate that the resulting equations 
for r and 7700,7710 are in agreement and produce results 
similar to those shown in Fig. [T] and O The expressions 
in Eq. (|48|) use the fact that the sums and integrals coin¬ 
cide in the continuum limit. To see that this is the case, 
note that the sums in ay have a factor of c more states 
for each interval of energy which cancel the factor of 1 /c 
multiplying each sum. 

Figure [5] shows the poo (black solid line) and 7710 /10 
(blue dotted line) from Eq. (H51) as a function of the ratio 
of y- to x-length scale at the threshold energy E = Ei^i. 
The 7710 stays positive for the values shown but the 7700 
changes sign for c ~ 4.89. At this value of c, the effect 
of the infinite closed channels disappears for any finite 
value of the scattering length. In other words, the scat¬ 
tering behaves as if there were no confining potential at 
all; a similar effect also exists in the harmonic oscillator 
confining potential which occurs at specific values of en¬ 
ergies away from the channel thresholds M- Over the 
range shown in FiglH] at large values of c the threshold 
value of 7700 appears to be proportional to y/c while the 
7710 appears to be a constant plus a function proportional 
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to l/\/c. 

Figure [ 6 ] depicts the energy dependence of 7700 and 
7710/10 parameters for different values of c. As c becomes 
larger, the corresponding Hamiltonian changes from a 
quasi-one dimensional to a quasi-two dimensional one. 
Since the 7700 diverges just below each threshold, these 
curves exhibit many points of divergence corresponding 
to each threshold. However, an overall pattern can be 
seen in panels (b) and (c). This pattern is clearer in 
the plot of 7710/10 [see red dashed line in Figl 6 l[c)] which 
does not diverge at threshold. In FigsU] (a) and (b) the 
c parameter acquires the values c = 5 and c = 40 . In 
these panels is observed that as c increases the system 
behaves more like it is only confined in the a:-direction 
and therefore those thresholds [see the magenta arrows at 
e = 1 / 4 , 9 / 4 , 25 / 4 ,... if FiglBJb)] dominantly character¬ 
ize the energy dependence in the corresponding 7700 func¬ 
tions. Moreover, FiglBJb) demonstrates that the multiple 
divergences of the 7700 parameter occur around an enve¬ 
lope curve. This overall behavior can be captured by 
defining an averaged 7700 according to the relation 

{Voo) = 777- / Vooiv)diy ( 49 ) 

which will give a smooth curve in the limit that the en¬ 
ergy spacing is smaller than Sv. 


3. Off-center scattering for a square well confining potential 

The above results for the infinite square and infinite 
rectangular wells are based on the assumption that the 
scattering center, i.e. the short-range potential I4, is 
placed at the center of the confining potential. In this 
section, the confining potential is considered to be an in¬ 
finite square well where the scattering center is located at 
a generic point, (xs,ys) = a±{S:s,ys) whereby the trans¬ 
verse wave function |$(xs, t/s)! is no longer only 0 or 2 /a_L 
as in the previous cases with {xs,ys) = ( 0 , 0 ). Note that 
similar systems have been investigated where the inter¬ 
acting particles are confined in separate harmonic waveg¬ 
uide geometries yielding in this manner confinement- 
induced interlayer molecules [s^. 

In the following, the discussion restricted on s—wave 
scattering which couples on even z-parity states of the 
confining square well and/or p-wave collisions which cou¬ 
ple on odd z-parity states of the corresponding confining 
potential. 

In this case the eigenenergies are defined as follows: 

h^TT^ 2 2\ 

Ena;ny = 2^0^ 

but now the rix and riy can independently be 
1 , 2 , 3 . The one dimensional wave function is = 

a±_ sin(na; 7 r[i — 1 / 2 ]). The transverse eigenfunction 
of the corresponding confining Hamiltonian at the 


scattering center behave as 
2 - 

^a{Xs,ys) = — sm{kx[xs-l/2])sm{kx[ys-i/2]) ( 51 ) 

a_L 


where kx = tLx'k and ky = riy-K. The only nonzero terms 
when the scattering center is at the origin is from the rix 
and riy which are the odd integers. 

For the case of s-wave scattering which couples on even 
parity states of the square well and for the case of p- 
wave scattering which couples on odd z-parity states the 
prescription given in the appendix 1 C 21 defining in this 
manner the corresponding 77-parameters. Namely, the 
7700 and 7710 functions fulfill the following relations: 


Voo 


—4 ^ [ Rnx,ny{XsTys) 


^ — ckq —eAz 


“n>e 


-4 


/*oo g —27r\/Q(^ —eAz ^ 






v'« 0 - 

ada 


Vw = 


with 


_ _ 

-487r^ ^ Rn^,ny{,is,ys)\Jal - e 

al>t 

poo 

+48 tt^ - e ( 52 ) 

J \/e 2 


Rn^,ny{Xs,y8) = -^\^a{x,y)\^ 

= sin^(fca:[i — 1 / 2 ]) siv?(ky[y — l/ 2 ]j[ 53 ) 


where ag = [nxj^'f' -|- (r7j^/2)^. In the integrals, the 
double integrals duxduy have been converted to cylindri¬ 
cal coordinates and the integral over angle from 0 to 'k 12 
has been carried out. It may appear that these expres¬ 
sions do not satisfy the condition that the sum and the 
respective integral should have the same form in order for 
their divergences to cancel and leave a finite difference. 
To see that the condition still holds, note that the sum is 
over 4 times as many terms as the centered square well, 
but the average value of Rnx,nyixs,ys) is 1 / 4 . Numer¬ 
ically tests have confirmed that the resulting equations 
for F and 7700,7710 exhibit excellent agreement, and they 
depict features similar to those shown in Figs. [T] and [2] 
Figure [ 7 ] depicts the 7700 function for different positions 
of the scattering center: {xs,ys)- As in Fig. HJ the quan¬ 
tity 7700 is singular as the energy approaches thresholds 
from below. In addition. Fig. [ 7 ] demonstrates that there 
are many more thresholds when the scattering center is 
shifted from the origin. This is due to the possibility of 
scattering into states that have a node at the origin. In 
discussing the position of thresholds, we will only give 
the values for rix < ny. The black solid line in panels 
(a) and (b) refers to the position (0,0) and is the same 
result as plotted in Fig. 01 At this position, there are 
only singularities at the values of rix , ny equal to odd in¬ 
tegers. The thresholds in Fig.[ 7 ]are 1,3 at e = 2 . 5 , 3,3 at 
e = 4 . 5 , 1,5 at e = 6 . 5 , and 3,5 at e = 8 . 5 . In Figl 7 ][b) 
the 7700 parameter is shown for a scattering center being 
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Figure 7. (color online) The parameter rjoo for infinite square 
well confinement is shown as a function of the scaled energy 
e = i5/(4ft^7r^/[2/ra5_]) with the scattering center at the scaled 
position {xs,ys)- Panel (a) (0.0, 0.0) is the black solid line, 
(0.1, 0.1) is the red dot-dash line, and (0.2, 0.2) is the green 
dot-dot-dot-dash line. Panel (b) (0.0, 0.0) is the black solid 
line, (0.1,0.0) is the red dotted line, (0.2,0.0) is the green 
dashed line. The solid line is the same result as in Fig'Hl 


placed along the x-axis, namely at (0.1,0.0) (red dotted 
line) and (0.2,0.0) (green dashed line) can have addi¬ 
tional singularities than in the case of a scatterer placed 
at the origin [see black solid line in FigITKb)]. Specifically, 
the additional structure emerges when one of nx,ny is an 
even integer. The additional singularities over those for 
the scatterer at the center are for 1, 2 at e = 1.25, for 2, 3 
at e = 3.25, for 1,4 at e = 4.25, for 3,4 at e = 6.25, for 2, 5 
at e = 7.25, and for 1, 6 at e = 9.25. In FigITKa) the ?7oo 
parameter is shown for a scattering center being placed 
along the diagonal of the x — y plane, namely at (0.1, 0.1) 
(red dot-dash line) and (0.2,0.2) (green dot-dot-dot-dash 
line) where both curves possess additional threshold sin¬ 
gularities in comparison with the black solid line and the 
curves in Fig. [7Kb). Specifically, the additional structure 
emerges when both nx,ny are even integers. The addi¬ 
tional threshold singularities are for 2, 2 at e = 2, for 2,4 
at e = 5, for 4,4 at e = 8 and for 2, 6 at e = 10. 

The coefficient of the diverging term is proportional to 
Rn^,ny and, thus, has different effect depending on the 
value of the transverse function at the scattering center 
that corresponds to that threshold. For example, the di¬ 
vergence at e = 2 (corresponding to 2,2) is clearly visible 
when the scattering center is at (0.2,0, 2) but is not vis¬ 
ible when the center is at (0.1, 0.1). This is because the 
latter position is not far shifted from the origin so that it 
is still near enough to the node of the transverse function 
to be suppressed. The divergence at e = 5 (correspond¬ 
ing to 2,4) is visible in both because the larger quantum 
number increases the amplitude of the transverse func¬ 
tion at the point (0.1,0.1). Another interesting feature is 
that the divergence at e = 4.25 corresponding to 1,4 has 



e 


Figure 8. (color online) The total absorption probability 
(in arbitrary units) for a harmonic oscillator confining 
potential as a function of the scaled energy, e = E/{fko±). 
The lines correspond in calculations using Eq. (|40|) with 
{i,m) = (0,0). The symbols refer to the ab initio quantum 
calculations using coupled angular momenta for the energy 
range 0.9 < e < 5.5. In panel (a) the s-wave scattering lengths 
take the values of —a±/4: (black solid line) and —a^l2 (red 
dotted line). In panel (b) the s-wave scattering lengths refer 
to values of ax/4 (black solid line), ax/2 (red dotted line), 
and 0.68475ax (green dashed line). 

a large strength partly because it is near the divergence 
at e = 4.5 corresponding to 3,3. 

There are similar types of features in the plots of rjiQ. 
However, because 7710 does not diverge at thresholds, 
there are discontinuities in slope at the newly allowed 
thresholds. 


B. Photoabsorption cross section in a harmonic 
and square well confining geometries 

This subsection investigates the impact of the transver- 
sally harmonic or square well confining potentials on the 
photoabsorption cross sections. In Fig. [SI the absorp¬ 
tion probability is plotted as a function of the scaled en¬ 
ergy e = E/fiu}± for the harmonic oscillator confining 
potential. The different plots are for different scattering 
lengths for {£,m) = (0,0). This figure has both lines 
and symbols. The lines are from Eq. (l40l) while the sym¬ 
bols are from an ab initio quantum calculation that used 
a grid in r and angular momentum to compute the to¬ 
tal outgoing flux. The only adjustment between the two 
types of calculations was an overall scale size. The good 
agreement between the two types of calculations shows 
that Eq. (gOl) is an excellent approximation. 

Two interesting features are worth noting. The most 
obvious feature is that the total probability goes to 0 
at every threshold. The probability goes to zero as the 
threshold is approached from below and from above. As 
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Figure 9. (color online) The partial absorption probability, 
1^ (in arbitrary units) for a harmonic oscillator confining 
potential in arbitrary units as a function of the scaled energy, 
t = E/{fkjj±). The lines are calculations using Eq. (1401) with 
{l,m) — (0,0). Panels (a) and (b) correspond to s-wave scat¬ 
tering lengths of —ax/4 and Ox/4, respectively. The black 
solid line is for absorption into the 1st open channel, red dot¬ 
ted line is 2nd open channel, green dashed line is 3rd open 
channel, and blue dot-dashed line is 4th open channel. 


Figure 10. (color online) The total absorption probability 
\D~ f (in arbitrary units) for an infinite square well confining 
potential in arbitrary units as a function of the scaled energy 
e = i?/(4fi^7r^/[2/rax]) with {£,m) — (0,0). The lines refer to 
calculations based on Eq. (001). In panel (a) the s-wave scat¬ 
tering lengths acquire the values —ax/16 (black solid line), 
—ax/8 (red dotted line) and — ax/4 (green dashed line). In 
panel (b) the s-wave scattering lengths have the values ax/16 
(black solid line), ax/8 (red dotted line) and ax/4 (green 
dashed line). 


the energy approaches a channel threshold from below, 
the parameter 7700 —>■ 00 which means that f 0 

since the ryoo is in the denominator of \D~\'^. As the en¬ 
ergy approaches a channel threshold from above, Nqq —>■ 
00 because the corresponding channel momentum van¬ 
ishes, namely fc /3 —>■ 0. In the expression for \D~\'^, there 
is an Nqq in the numerator and an Nqq in the denomina¬ 
tor yielding \D~\^ 0. The total absorption probability 

is proportional to ^/lE — Ea\ at each channel threshold 
Ea ■ This behavior is in strong contrast to the unconfined 
absorption probability where there is a discontinuity in 
the slope at each channel threshold, but the total ab¬ 
sorption does not go to zero at each threshold. Also, in 
the case of confinement, the photoabsorption probabili¬ 
ties possess discontinuous slopes near the channel thresh¬ 
olds. 

Another important feature is that a negative scatter¬ 
ing length causes a resonance-like structure just below 
each channel threshold [see FigjSKa)]. On the other hand 
in FigjSKb) no resonant structure emerges below each 
threshold if the s-wave scattering length is positive. This 
behavior can be understood by the fact that for nega¬ 
tive s-wave scattering lengths the term 1 -|- as{E)r]QQ{E) 
in the denominator of the total absorption probability 
vanishes yielding in this manner a maximum in the total 
probability. This effect clearly originates from the closed 
channel physics. In addition, in FigjSKa) we observe that 
for smaller (in magnitude) values of scattering length, 
the resonance occurs closer to threshold and is narrower 
in energy. This is to be expected because the smaller 
(in magnitude) scattering length corresponds to a state 


that is more weakly bound and, hence, a smaller overlap 
with the scattering center leading to a longer lifetime and 
narrower resonant features. 

Figure O depicts the partial absorption probability 
for a negative [see Figl^Ka)] and a positive s-wave 
scattering length [see Fig[9Kb)] case. The partial cross 
section is proportional to |Wa^ooP which is proportional 
to l/ka since the <I>o,(0,0) has no dependence on the a 
channel quantum number. The channel momentum ka 
decreases with increasing a which means the partial cross 
section increases with increasing a. Also, the partial ab¬ 
sorption cross section exhibits differences between chan¬ 
nels just above a channel threshold as it is demonstrated 
in Fig. E i e. in panels (a) and (b) for the energy interval 
3 < e < 5 the \D~\‘^ of the second channel (red dotted 
line) is larger than the partial photoabsorption in the 
first channel. In this case, the channel that just opened 
will have ka ^ 0 and will therefore have nearly all of the 
outgoing probability. 

Figure [TOI demonstrates the cases of total absorption 
probabilities for the infinite square well case for negative 
[see Fig. unia)] and positive [see Fig. [TUT b)] s-wave scat¬ 
tering lengths. The black solid lines indicate the s-wave 
scattering length with magnitude ax/16, the red dotted 
lines refer to magnitude a±/8 and the green dashed lines 
denote magnitude a±/4:. The general features are the 
same as for Fig. [S] Unlike the Fig. [ 8 | case, there are 
degenerate channel thresholds for the square well confin¬ 
ing potentials. For example, the first excited threshold 
at e = 2.5 is doubly degenerate: nx^Uy of 0,1 and 1,0. 
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Figure 11. (color online) The total absorption probability 
(in arbitrary units) for a harmonic oscillator confin¬ 
ing potential is plotted as a function of the scaled energy, 
e = In panel (a) and (b) the lines are calcula¬ 

tions using Eq. (gOl). The black solid line is the result when 
Voo = 0. The red dotted and green dashed lines correspond to 
as{E)/a± = 0.02 and —0.02 respectively. The blue dash-dot 
and purple dot-dot-dash-dot lines are as{E)/a± = 0.04 and 
—0.04 respectively. Panel (b) consists of a magnification of 
panel (a) around the threshold energy e = 3. 


Also, not all thresholds are equally spaced, in particu¬ 
lar for e > 10 which leads to features that change with 
threshold. 


C. Negative-ion photodetachment in magnetic 
fields 

The result of photodetachment of an electron from a 
negative ion in a magnetic field was treated in Ref. @ 
and is similar to the example of photodetachment in 
an isotropic, two-dimensional harmonic oscillator poten¬ 
tial. However, in Ref. only the effect of open chan¬ 
nels is taken into account. This is equivalent to set¬ 
ting r]£m^) = 0 in Eq. (l40l) above. Thus, the results 
in Ref. are an approximation to the full treatment 
which includes the effects of the closed channels. Since 
the formalism of Ref. @ has been successfully applied in 
many circumstances, it is worth investigating the regimes 
where this approximation {rjirn = 0) is adequate. Figure 
[TT] depicts the total photoabsorption probability (in ar¬ 
bitrary units) for different s-wave scattering lengths that 


are a small fraction of a_L which is the typical case for 
photodetachment in laboratory strength magnetic fields. 
The black solid line is the result from the approxima¬ 
tion of Ref. @ . For small positive scattering length (see 
red dotted and blue dash-dot line in Fig fTTT a)) , the pho¬ 
toabsorption cross section vanishes below and above each 
channel threshold whereas the black solid line vanishes 
only above every the threshold. Note that the cross sec¬ 
tion vanishes around every threshold also for negative 
scattering lengths [see green dashed and purple dot-dot 
dashed lines in Fig fTlT a)]. Moreover, the case of negative 
scattering lengths exhibits another qualitative difference 
from the approximation of Ref. @ . More specifically, this 
can be seen in Fig fTTT b) where the total photoabsorption 
probability for negative scattering length (green dashed 
and purple dot-dot-dashed lines) possesses a resonance 
barely below the channel threshold. This non-trivial res¬ 
onant feature is absent from the black solid line manifest¬ 
ing in this manner the importance of the closed channel 
physics. In addition, FiglTT](a) and (b) shows that away 
from the channel thresholds our results are practically the 
same as in the approximation of Ref. @ , [black solid line 
in Fig in]- Note that there is one energy between each 
threshold where all of the results are the same. These 
are the energies where 7700 = 0 . 

Evidently, the role of the closed channels becomes im¬ 
portant only around each threshold due to small scatter¬ 
ing lengths. This explains why the result of Ref. Q works 
well for photodetachment in a magnetic field. In this sit¬ 
uation, the scattering length is of the order of a few Bohr 
radii while the confinement length is a_L = ao\/ Ba.u./B 
where oq is the Bohr radius and Ba.u. is the atomic unit 
of magnetic field. Thus, even in a 1 T magnetic field, the 
scattering length will typically be less than a_L/100. How¬ 
ever, we expect that experiments would show differences 
from the treatment of Ref. @ in stronger magnetic fields 
or for a negative ion whose photodetachment probes a 
final state having a much larger electron-atom scattering 
length. 

The photodetachment for higher partial waves may 
possess shape resonances which allows for a larger ef¬ 
fect from the closed channel collective state described 
by rjirniE). Therefore in the following, consider a pho¬ 
todetachment process into i = l,m = 1 for a harmonic 
oscillator confinement. The confining potential does not 
have as strong an effect on the £ = l,m = 0 case because 
the outgoing flux is mostly along the z-axis compared to 
m = 1 where the outgoing flux is mainly perpendicular 
to the z-axis. The energy dependent scattering volume 
is given by the relation iaii5i{E)/[ka±Y where a reso¬ 
nance form for the phase shift, (5i is chosen according to 
the relation: 


tan^i = —k^ 



e — Cr _ 


(54) 


where e = E/{huj±) is the scaled energy with E = 
h^k"^/2^, Cr = Erl is the scaled energy of the reso¬ 
nance with Er = h^kr/2fi, 7 is a parameter proportional 
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Figure 12. (color online) The I = \ phase shifts used in the 
calculations of Fig. ll3l The black solid line is a non-resonance 
example and the red dotted line is the resonance case. The 
energy of the resonance is tr = 7.8. 


to the scaled energy width of the resonance, and fee is 
the scaled scattering volume. Figure [T^] shows the (5i 
phase shifts for a resonance case (red dotted line) with 
fsc = —(1/8)^, 7 = 1.4, and e^. = 7.8 and for a non¬ 
resonance case (black solid line) fsc = —(1/8)^ and 7 = 0. 
For the energy dependent smooth dipole, Ds{E), the fol¬ 
lowing form is chosen 


Ds{E) 




ds + {dr/a]_) 


7 /^ 

e — Cr 


(55) 


where dg represents the amplitude for photoabsorption 
directly to the continuum and dj. represents the ampli¬ 
tude for photoabsorption into the resonance. Figure [T51 
shows the results from 4 different calculations. In all cal¬ 
culations, the green dotted line is the photoabsorption 
that would occur if there were no conhning potential, the 
black solid line corresponds to the full calculation with 
the confining potential, and the red dashed line is the ap¬ 
proximation in Ref. Q where the effect of the closed chan¬ 
nels is neglected which is equivalent to setting ryu = 0. 
As expected, the resonance features from the confinement 
calculations give an overall shape that tracks the non¬ 
confinement calculations, but with sharp features near 
the thresholds. The topmost calculations are when there 
is no resonance 7 = 0. The “Lorentz” calculation has 
ds = 0; it does not look like the standard Lorentzian 
shape resonance due to the overall increase proportional 
to The “Window” calculation has dr =0; as with 
the Lorentz calculation, this does not look like the sym¬ 
metrical window resonance due to the overall increase. 
The “Fano” calculation has dr/dg = 60 which gives a 0 
in the dipole matrix element at e ~ 9.1 which is at a 
somewhat higher energy than the resonance position at 
7.8. The full and approximate confinement calculations 
differ most strongly where tan di is large which is near the 
resonance and an energy range near the thresholds. For 
the three cases with the resonance at er = 7.8, the ab¬ 


sorption cross section shows sharp resonances below each 
threshold only for e < tr because the scattering volume 
is negative in this range. 



Figure 13. (color online) The ^ = 1 photoabsorption 
(in arbitrary units) for different cases. In all calculations, the 
green dotted line is for photoabsorption with no confining po¬ 
tential, the red dashed line is the approximation that includes 
the confining potential but does not include the effect of closed 
channels (equivalent to setting 711 = 0 ), and the black solid 
line is the full calculation including the confining potential. 
The “Smooth” calculation does not have a resonance. The 
“Lorentz”, “Window”, and “Fano” calculations only differ in 
the energy dependence of the dipole matrix element when no 
field is applied. 


VII. SUMMARY AND CONCLUSIONS 
A. Summary of Method 

It is evident that there are many steps in the deriva¬ 
tion and application of the Schwinger variational ap¬ 
proach to collisions involving a short range spherically 
symmetric potential in the presence of various types of 
confining geometries. Therefore, in this section, the im¬ 
portant steps are pointed out for applying the method. 
First, Eq. (1331) gives the expression for the AT-matrix and 
Eq. (1361) gives the dipole matrix element in terms of the 
scattering length or volume and the energy dependent pa¬ 
rameters Uajm and rjim- The conversion from spherical 
coordinates to the confinement geometry is encapsulated 
in the Ua,im in Eq. (1341) which depends on the conhne- 
ment wave functions <I>Q,(a;,y) at the scatterer. 
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In this paper, the expressions for the r]im were obtained 
from a detailed examination of the confinement Green’s 
function, Gc, and various integral expressions for the free 
Green’s function, G /. There is a conceptual shortcut that 
can be used to obtain the ijim- the Gc becomes the G/ 
in the limit that the confinement length scale, a_L, goes 
to infinity. This is the reason that the r]im have the form 
of a sum (from the Gc) minus an integral that has the 
same functional form (from the G/ as formulated from 
Gc in the limit a_L —>■ oo). Similar to Eq. (l34l) . define a 
transformation for the closed channels 




,im 


V67rKQ,a_L$a(ps) 

.(a±$o) 


i = m = 0 

e=l, m = 0 (56) 

£ = 1, m = ±1 


where the operators are defined below Eq. (051) , Ps is the 
position of the scatterer, and E = —h‘^K^/{2pa\) + Ea- 
Since is proportional to l/aj_, the U only depends 
on the form of the confining potential, the scaled energy, 
and the Pstaj_. In this paper, we restricted the confining 
potential to be independent of z and symmetric in xy. 
For these cases, the following relation between y and U 
hold 

OO 

mm = ^ 

Oi'xy.o 

- / (57) 

J (y.>oco 

where the limit Az —>■ is understood. The integral 

over da is the result of the quantum numbers for G/ 
giving an infinitesimal spacing of Kq, in the limit a± —>■ 
oo. For the case of Sec. I VIA 31 which does not satisfy 
the restrictions on the confining potential, the \Ua,im'Y is 
replaced with the average value as a_L —> oo. 


B. Conclusions 

In the preceding sections, the development of a 
Schwinger variational framework is presented as a treat¬ 
ment of scattering in a confined geometry. The current 
formalism is non-perturbative, permitting the treatment 
of a wide class of Hamiltonian systems which possess two 
potentials. The particular type of Hamiltonians stud¬ 
ied in this work consists of systems where the inter¬ 
action at short distances is spherically symmetric and, 
at large distances, a confining potential is considered 
which bounds the motion of the particle in the degrees 
of freedom perpendicular to the direction of its propaga¬ 
tion. The restriction in the current formalism amounts 
to the fact that the considered potentials dominate at 
different scales. For the sake of simplicity, we restricted 
the confining potential to have no dependence in the z- 
coordinate and symmetry in xy. However, this theoreti¬ 
cal method could treat general confining potentials where 
length scale separation holds. 


The theoretical formalism presented in this paper al¬ 
lows a self-consistent treatment of scattering in a confined 
geometry where the results are manifestly convergent at 
each step of the derivation. This formalism has been 
used to derive results obtained in previous studies (e.g. 
s- and p-wave scattering within a harmonic confining po¬ 
tential and s-wave scattering in a square well confining 
potential) and to quickly derive results in novel geome¬ 
tries (e.g. s-wave and p-wave scattering in off-center con¬ 
fining potentials and scattering in rectangular confining 
potentials). Results for the scattering parameters in dif¬ 
ferent geometries were presented and main qualitative 
features were explained. The formalism was also used 
to derive a treatment of half-scattering problems (e.g. 
photo-detachment) in confining potentials. For example, 
the photo-detachment results in a harmonic confinement 
was compared to a previous LFT treatment that only ac¬ 
counted for open channels to assess the role of the closed 
channels which are included in this treatment: the ef¬ 
fect from the closed channels are most important near 
thresholds and become more important as the scatter¬ 
ing length becomes a sizable fraction of the confinement 
length scale. 

An interesting possibility for the formalism in this pa¬ 
per is to aid in the understanding of the implementa¬ 
tion and possible limitations of the LFT. The basic steps 
in the current formalism were cast in a form strongly 
reminiscent of the LFT. However, these two frameworks 
possess a conceptual difference which is mainly focused 
on the physics of the closed channels. More specifically, 
in the LFT approach the AT-matrix formulas include di¬ 
verging sums and therefore regularization schemes are 
employed to remove such singularities. In the Schwinger 
variational approach, such divergences do not emerge, 
which eliminates the need for tricky regularization tech¬ 
niques that are often used in the LFT theory. In particu¬ 
lar, we observe that in the physical LC-matrix of the LFT 
the information of the free space G/ is absent and it has 
been indirectly incorporated via regularization schemes 

[IlHlIl- 

In view of the rigorousness of the Schwinger variational 
approach, we expect that it can be equally applied to 
other physical systems where the LFT method has been 
used. For example, this method can be used to derive 
the scattering of an electron in a Rydberg state from 
a neutral perturber and, for £ = 0, the same result as 
Eq. (14) of Ref. [ 2 ^ is obtained. Note that the frame 
transformation ideas have been used in many different 
circumstances. A short list of examples could include the 
jj — LS transformation, molecular rotational transforma¬ 
tions, molecular vibrational transformations, local frame 
transformations involving external electric and/or mag¬ 
netic fields. In all cases, wave functions in one represen¬ 
tation are projected on those in another representation 
at a surface (or indirectly in a region) where both rep¬ 
resentations are expected to be accurate. However, the 
level of error involved in such a procedure is not clear. 
It is conceivable that the method described above might 
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enable a more systematic derivation of the frame trans¬ 
formation so that the level of expected error would be 
clearer. 
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Appendix A: Matrix elements with E, IV'.} 

The volume integral in Eq. (I25|) can be simplified by 
employing the relation 14 |^s) = [E — Then, 

assuming length scale separation, the confining potential 
14 is practically zero in the range of 14 where the follow¬ 
ing relation is valid {E — He) |'0c) ~ {E — Hf) j'i/i’e) ~ 0. 
Thus, integrating by parts twice Eq. (I2S1), the terms of 
[E — Hf) \'tl)e) vanish and only the surface ones survive. 

In general, the following relation is fulfilled: 

{'4ic\VsHs) = ^ [ n-{^*{r)[Vi^s{r)]-[V^*{r)]iPsir)}da 

Ja(T) 

(Al) 

where cr(r) indicates the surface containing the volume 
T and h is the unit vector which is outward normal to 
the surface. If the region r is a sphere of radius r, this 
matrix element can be written as 


term. However, at small distances only few terms are 
needed to be matched. 

The coefficient of the jo{kr)Yoo{^) term is simply 
times the value of '0c (r) at the origin giving 

0c(r) - i^ciO)'^joikr)Yoo{n) (A3) 

After substituting this expression into Eq. (IA2I) the sur¬ 
face integral yields the following general form: 


(0cl^s|0s) = 


2/r 


0:(r)^-^0,(r)JdD 

(A2) 

where dD = sin 6d6d(j) is the solid angle differential ele¬ 
ment. 

This subsection recast the integrals which involve 
14 10s) into surface integrals when length scale separa¬ 
tion between 14 and 14 is a good approximation. The 
following subsections focus on special cases where |0s) is 
an eigenstate of angular momentum i and its projection, 
providing explicit analytical equations for this situation. 


(0c|t4|0s) = y ^0c(O)y^ ^-itando^ , (A4) 

where the expressions have been used for the 0c (r) = 
0c?2(r) [see Eq. (O] and 0s(r) = 0s,oo(i’) Eq- (|A4l) reads 

= Ki0,0)sl^ (-^tando) (A5) 


2. Analytic expressions of (0c|Vs|0s) for £ = 1, m = 0 

This subsection focuses on the case where |0s) = 
10s,^m) possesses an £ = 1, m = 0 character. In a 
similar way as was discussed in the previous subsection 
the 0c (r) functions are matched to the superposition 
of spherical Bessel functions times spherical harmonics, 
ji{kr)Yim(yt)- In this particular case the coefficient of 
the ji(A:r)Yio term is simply jk times the value of 
90c(r)/c}z at the origin. 


0 c (r) ~ 


50c (r) 


dz 


ji{kr)Yio{n) 


|r|=0 


Substituting this expression into Eq. 
integral reads: 


(A6) 
the surface 


(0c|Es|00 = 


I nhP 50c (r) 


d 


dz 


|r|=0 


IQtt f 1 


(A7) 

where using the expressions for the ipciy) = 0c(a(r) from 
Eq. ([5]) and 0s (r) = 0s,io(r) Eq. (IA7I) reads 


(0(“)|t>s|0s.io) = $:(o,o)y^ 


67rfco 

k^ 


—-i-tandi 1 (A8) 

TT 


1. Analytic expressions of (0c|t4|0s) for f = 0, m = 0 

Assume that in this case study the states |0s) = 
10s.fm) possess I = Q, m = 0 character. At distances 
where the 14 potential is considered small, the 0c(r) 
to a good approximation can be written as a linear 
combination of spherical Bessel functions times spheri¬ 
cal harmonics, ji(kr)Yimi^)- Note that the wave num¬ 
ber k corresponds to the total energy and G denotes the 
{6, 0)-angles. The coefficients of the linear combination 
can be specified by expanding in Taylor series the 0c (r) 
and ji{kr)Yimi^) functions and matching them term by 


3. Analytic expressions of (0c|t4|0s) for 
f = 1, m = ±1 


This subsection considers the case where |0s) = \tps,im) 
possesses £ = 1, m = ±1 character. Therefore in this 
case the coefficient of the ji{kr)Yi±i(n) term is simply 
—\J2Yk jk times the value of (50c(r)/5a:=F*50c(r)/5t/)/2 
at the origin yielding the following relation: 


0c(r) ~ -(5±0c) 



ji(fcr)Yi±i(D) 


(A9) 


|r|=0 
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where (9±'!/'c) = /dx ^ idipci^)/dy)/2. 

Substituting these expressions into the surface integral 
gives the simple result: 


0 

(AlO) 

where using the expressions for the ipc(j) = from 

Eq. (O and V's(i‘) = '0s,i±i(r) Eq. (lAlOl) reads 


= (B3) 

Similarly, for the case of£=l,m = ±l the dominant 
contribution to the Ds matrix elements emerges from the 
term -pi ■ p2AGpp{0) 

Ds,i±i = - {iJs,i±i\Vs\pie^^'^^) AGppiO) {p2e^^'^^\Vs\'4’s,i±i) 

= (B4) 


{-IpclVsl-lps) = 


/ nti^ 


id±tpc)* 


{i^^j:i\Vs\^s,i±i) 


-(5±$a)* 



-tan 

TT 


<5i 

(All) 


Appendix B: Evaluation of the Ds matrix elements 

At this point, all of the terms in the expression for 
the AT-matrix, Eq. dSD have a simple analytic expres¬ 
sion except for the term involving the matrix element of 
VgAGVs, see Eq. The most straightforward method 
for evaluating this matrix element is to individually com¬ 
pute the matrix elements of VsGfVg and VsGcVs and sub¬ 
tract. However, this method is fraught with difficulties 
because each of these matrix elements involve integrands 
that diverge. 

At the symmetry point of Vc, the function AG only has 
even powers of {zi — Z 2 ) and powers of pi ■ p 2 through 
second order if the confining potential is a power series 
in x,y,z. The vector pi = xxi -I- yyi. Thus, the matrix 
element in Eq. (1241) only has terms from regular functions 
at least to order r^. If the state Itps) possesses £ = 0 or 
£ = 1, m character, then only terms in AG proportional 
to {zi — 2 : 2 )°, {zi — 22 )^, and (pi — ^ 2 )^ will contribute 
to the integral: 


AG(ri, r 2 ) ~ AG(0) - AGzz{0)ziZ2 - AGpp(0)pi • p 2 , 

(Bl) 

where AG(0), AGzz and AGpp are coefficients that will 
be determined below. Note that terms proportional to 
xj, yj, and 2 ? are not included since they only affect 
matrix elements for d-wave scattering or give small cor¬ 
rections to the s-wave scattering. Recall that AG is a real 
function since its constituents (i.e. G/ and Gc) are the 
principal value Green’s functions. Using this expansion, 
the dominant contribution to the Dg matrix element for 
the £ = 0, m = 0 case is the AG(0) term 

Ds,oo = {'iPs,oo\Vs\z°i) AG{0) {z°\Vs\i;s,oo) 

= AG(0) — ^ f - - tan do) (B 2 ) 

p k \ TT J 

For the case of £ = 1, m = 0 the dominant con¬ 
tribution to the Dg matrix elements is from the term 
-ZiZ2AGzz{,Q) 

Ds,W = - {ll}g,w\Vs\zi) AGzz{^) {Z2\Vs\ll}g,w) 


Appendix C: r/tmiE) functions for given confining 
geometries 


In the expressions for the Green’s functions, we only 
need terms through pi - P2- The free space Green’s func¬ 
tion, G/ is Taylor expanded up to terms proportional to 
pi ■ p2 which give 


cos(fcri2) 

|ri2| 


cos(fcAz) 

|A5| 


Pi ■ p 2 


cos{kAz) 

|A5|3 


k sin(fcAz) 

Ai2 


(Gl) 

where the scaled variables are k = ka± and A 2 = Az/a_L. 
The strategy below will be to write this expression for G/ 
as an integral with a uniformly converging integrand. We 
will show that this integral will be the continuum form 
of the sums that arise in Gc- 

In this subsection, the functions rje-miE) are derived 
for two types of confining potentials, i.e. harmonic and 
infinitely square well confining potential. In addition, 
we compare the approximation Eq. dsm to the full ex¬ 
pression for r(£') as a function of Az, Eq. (IHTl) as an 
illustration of the accuracy of the rjimiE). 


1. Harmonic Oscillator Confining Potential 


In this section, the motion of the particle in the 
transversal degrees of freedom is bounded by a harmonic 
oscillator which possesses the form Vc(p) = ^u!±p'^ where 
p denotes the mass of the particle, uj± corresponds to the 
frequency of the harmonic confinement corresponding to 
the length scale aj_ = i/hJJuJZ- Note that the variable p 
represents the polar coordinate, namely p = \/x"^ + 
Under these considerations, the motion of the photo¬ 
electron separates in the 2 , p degrees of freedom. The 
solutions of the corresponding Hamiltonian, namely He, 
possess the form given in Eq. ([5]) for even or odd par¬ 
ity in ^-direction. Accordingly, the eigenenergies of the 
Hamiltonian He fulfill the relation E = huj±{2n -I- |m| -|- 
1) -|- 1 ^ 1 /(2/r), where n = 0,1,.... Evaluating the 

harmonic oscillator eigensolutions at the {p, <(>)—>■ 0 the 
following relations are obtained 


f$n,m(0) = 


a±^/n 


for m = 0 
for m = ±1 


(G2) 
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The scaled energy is defined by the relation e = 
E/{huj±). The scaled momenta can be written 
as fc = knm = \/2e — 2(2n + |m| + 1), and 

knm = \/2{2n + \m\ + 1) — 2e. Substituting Eq. (IC2I) 
in Eq. (I5T]) with the approximation of Eq. (ICll) . the fol¬ 
lowing relation is obtained 




+ 2 


n—0 n—no ,0 + 1 

rno,o sin(fc„oA.z) 


-dn — 2 




^nO 


•2—1/2 ^nO -tno.o 

- + 1) sin(fc„i A^) 

_p,.p24^-f- 


^nO 


dn 


n—0 


where the limit Az —>■ 0^ is understood for all expres¬ 
sions. Notice that all of the terms have the form of a 
sum minus a corresponding integral. This feature en¬ 
sures that the are well defined for any finite A5 and 
thus they give unambiguous, finite results in the limit 
A5-^ 0+. 

To summarize, the expressions needed for the calcu¬ 
lation of the A"-matrix are convergent series when using 
a finite Az in the expression for the Green’s function. 
When the scattering potential has a short range, ana¬ 
lytic expressions can be obtained in the limit Az —>■ 0+. 


2. Infinite Square Well Confinement 


- - . (n-b 1) sin(fc„iAz) 

+ Pi ■ P24 / ^- -dn 


+ Pi • P24 


knl 

(n-b 


n—Uo 1 + 1 


- Pi ■ P24 


(n -b l)e 


^nl 

— Knl A5 


-dn 


'rio,i 


^nl 


(C3) 


where no,m = (e— 1 — +|)/2 is the divide between open 
and closed channels; in the sums, the integer part should 
be used. Note that Pi ■ P 2 = pi ■ P 2 /This expres¬ 
sion has many terms but each can be identified with dif¬ 
ferent parts of the Eq. m- The first line is the con¬ 
tribution from the m = 0 open (first term) and closed 
(second term) channels from Gc- The second line is an 
integral expression that exactly equals —cos{kAz)/Az 
(see Eq. m) which is the contribution to m = 0 from 
Gf. The third (fifth) line is the contribution to m = ±1 
from the open (closed) channels in Gc- The fourth and 
sixth lines are an integral expression whose sum exactly 
equals the term in Eq. (ICll) multiplying the pi ■ P 2 - 
There are important features of Eq. (ESI) that affect 
how the scattering information is extracted. For exam¬ 
ple, the sums and integrals that correspond to the open 
channels all contain sin(fcQ,Az) and therefore they do not 
contribute to T as Az —>■ 0+ for Ap -C 1. Terms which re¬ 
fer to the closed channels do survive in this limit and they 
contribute to the r]im{E). Lines 4 and 6 contribute to the 
rji±i while the second term in each of lines 1 and 2 con¬ 
tribute to the 7700 and 7710 . The 7700 is the limit as all spa¬ 
tial coordinates of E go to 0. The 7710 = —3d^T/d(Az)^ 
in the limit all spatial coordinates go to 0. This gives 


g-K„oAz /-OO g-S„oAz 

7700 = 2 ^ ;-2 / —r- dn 


’/ 

'J Tin 


^10 


00 

= -6 ^ + 6 

-- I 1 n 


m±i —12 ^ 


n=no,0 + 1 ‘^ 0,0 

1 K-nl Az 


{n -b l)e 


^nl 


-12 


n=no,i + l 

(n-b l)e-”'‘i'^^ 


dn 




^nl 


(C4) 


In this section, the motion of the particle is bounded 
in the x- and 7/-direction by an infinitely square well. 
The square well possesses the following form V (x, y) = 0 
when |a;| < a_L/2 and \y\ < axl2- and V{x,y) = 00 else¬ 
where. As in the case of the harmonic oscillator at large 
distances the wave function of the electron separates in 
Cartesian coordinates, where the corresponding solutions 
of the Hamiltonian He possess the form in Eq. ([S]) for 
even or odd parity in the z-direction, respectively. In ad¬ 
dition, the energy spectrum of the He Hamiltonian obeys 
the following relation E = E^^^uy + k^k‘^^ ,^^/{2p) where 

^rixTiy = f^e 771 = 0 case, we are only 

interested in the values of the quantum numbers n that 
give non-zero functions at the origin; these correspond to 
the cases nx, ny = (2ax + 1), (2ay -b 1) where the ax, ay 
are independently 0,1, 2,.... For the tti = I case, we use 
the functions 




= -^[FnAx)Fny{y)±iEn,{y)Fny{x)] (C5) 


where Fn{x) = •\/2/a_L sin(n 7 r[i — 1/2]) with x = x/ax 
and nx,ny = (2ax + 2),{2ay -b 1 ) with the ax,ay in¬ 
dependently 0,1,2,.... This choice gives F!^^(0) ^ 0 
and Fn (0) 7 ^ 0. This particular linear combination in 
Eq. (IC5|) ensures that collisions with 777 = 1 affect the 
square well states whose x-dependent factor vanishes at 
the origin while the 7 /-dependent factor is nonzero, or vice 
versa. Evaluation of the square well eigensolutions at the 
origin gives 


2(-l)° 


|4>«,a„(0,0)= 

l(a±$± )| = 

kv ± aj,a„2|0 a_ 


— for 777 = 0 

c+o,y + i 2T7{a,, + l) 


a± 


for 777 = 1, 


where the functions (x, 7/) are defined in Eq. 


(C6) 


We define the 
E / / [2pa\]). 

can be written as 


scaled energy to be e = 

The scaled momenta 

k = 27ry^, kajn = 

2'K^e- [{ax + [|777| -b l]/2)2 -b {ay -b 1/2)^], and 

kam = + [ + 1 -b l]/2)^ -b {oLy -b 1/2)^ — e. 

Substituting Eq. (IC6I) in Eq. (I5T|) the following relation 
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is obtained: 


“ 0,0 . ,r . oo 

rt TP\ o sin(A:aoA2:) 

r(ri,r2 ,S) = - Stt > --^+87r > 

kaO 


q:=q:o,0 + 1 


+ Stt 


/ a, 

-1 


“‘’’0 sm{kcoAz) 2 Q 
-^- -a a — Sn 


- 1/2 


kct 


L 


f^aO 

oo —KqqAz 

^ t2 

—-- a a 


^aO 


- Pi ■ p232'K^ ^2 


3 ^ (ox + 1)^ sin(fcc«iA2) 


+ pi ■ P2S2 tv' 


L 


ka 

{ax + 1)^ sin(fcaiA2) 2 


+ pi ■ P232-7r® ^ 

q: = q:o 1 + 1 


1 ^al 

(a.+ 1)^6-"-^^' 


d a 


- pi ■ P2327r' 


p OO 

7 

./ Cio 1 


KqI 

\2 _ —ftry 1 Az 


(“2: + 1 ) e '^“1 ,2 


K.Q;1 


d^a (C7) 


where ao,m denotes the divide between open and closed 
channels and is given by (ax + [|ni| + l]/2)^ + (ay + l/2)^ = 
e; in the sums, the integer part of ax,ay should be used. 
Note that pi- P 2 = Pi' P 2 / ci\- This expression has many 
terms but each can be identified with different parts of 
the Eq. m- The first line is the contribution from the 
m = 0 open (first term) and closed (second term) chan¬ 
nels from Gc- The second line is an integral expression 
that exactly equals — cos(fcAz)/A 2 (see Eq. (ICip i which 
is the contribution to m = 0 from G/; the lower limit 
of integration is —1/2 for both and ay. The third 
(fifth) line is the contribution to to = ±1 from the open 
(closed) channels in Gc- The fourth and sixth lines are 
integrals whose sum exactly equals the term in Eq. m 
multiplying the pi - p 2 ] the lower limit of integration is — 1 
for ax and —1/2 for ay. Unlike the integral expressions 


for the harmonic oscillator, it is somewhat tricky to show 
that the integrals exactly equal the expressions from Gf. 
The simplest method we found was to change variables 
Sx = otx + {\'m\-\-\)/2 and Sy = ay-|-l/2 and then convert 
to cylindrical coordinates -I- and tan(0) = y/x 

and remember to only integrate from 0 < </ < 7r/2 since 
the Sx > 0 and Sy > 0+. 

To obtain the scattering parameters pirm we can re¬ 
peat the logic of the preceding section for the harmonic 
oscillator. This gives 


yoo 


Vio 




“0,0+1 


^aO 


— Stt 


/ 

«/ a 


00 —K,qqAz 

^ i2 

—r- d a 


00 p Q 

= —247r + 247r / 

_ I 1 d ar 


“0,0+1 


»7i±i = 96"^^ Y1 
“ 0,1 + 1 


(oix + 


—967r' 


poo 

I 

J OLq 1 


ACq-I 

( ^ I 1 \2 ^ —Koi Az 

(ax + 1 ) e 2 

- z - d a 

Kal 


~ —ftrvnAz i2 

Kaoe “ da 


(C8) 


where the limit Az —>• 0^ is understood for all expres¬ 
sions. Notice that all of the terms have the form of a 
sum minus a corresponding integral. This feature en¬ 
sures that the rjim are well defined for any finite Az and 
thus they give unambiguous, finite results in the limit 
Az —>■ 0. 

To summarize, the expressions needed for the calcu¬ 
lation of the A'-matrix are convergent series when using 
a finite Az in the expression for the Green’s function. 
When the scattering potential is short range, analytic 
expressions can be obtained in the limit Az —>■ 0+. 
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